Search for: All records

Abstract Eruptive mass loss in massive stars is known to occur, but the mechanism(s) are not yet well understood. One proposed physical explanation appeals to opacity-driven super-Eddington luminosities in stellar envelopes. Here, we present a 1D model for eruptive mass loss and implement this model in theMESAstellar evolution code. The model identifies regions in the star where the energy associated with a local super-Eddington luminosity exceeds the binding energy of the overlaying envelope. The material above such regions is ejected from the star. Stars with initial masses of 10−100M_⊙at solar and SMC metallicities are evolved through core helium burning, with and without this new eruptive mass-loss scheme. We find that eruptive mass loss of up to ∼10⁻²M_⊙yr⁻¹can be driven by this mechanism, and occurs in a vertical band on the H-R diagram between $3.5\lesssim \log \left(T_{\mathrm{eff}}/\mathrm{K}\right)\lesssim 4.0$. This predicted eruptive mass loss prevents stars of initial masses ≳20M_⊙from evolving to become red supergiants (RSGs), with the stars instead ending their lives as blue supergiants, and offers a possible explanation for the observed lack of RSGs in that mass regime. 

Abstract We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range toZ= 10⁻⁶toZ= 0.2, which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity rangeZ= 10⁻⁴toZ= 0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions betweenX= 0 andX= 0.1 includingX= 10⁻², 10⁻³, 10⁻⁴, 10⁻⁵, 10⁻⁶that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables inMESAand find that calibrated solar models agree broadly with previously published helioseismic and solar neutrino results. We find differences between using the new 1194 OPLIB opacity tables and the 126 OPAL opacity tables range from ≈20% to 80% across individual chemical mixtures, up to ≈8% and ≈15% at the bottom and top of the solar convection zone respectively, and ≈7% in the solar core. We also find differences between standard solar models using different opacity table sources that are on par with altering the initial abundance mixture. We conclude that this new, open-access set of OPLIB opacity tables does not solve the solar modeling problem, and suggest the investigation of physical mechanisms other than the atomic radiative opacity. 

Abstract We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The newauto_diffmodule implements automatic differentiation inMESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite-difference approximations. We significantly enhance the treatment of the growth and decay of convection inMESAwith a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron-degenerate ignition events. We strengthenMESA’s implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars inMESA, we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator-split nuclear burning mode. We close by discussing major updates toMESA’s software infrastructure that enhance source code development and community engagement. 

Abstract We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729M_⊙carbon–oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of²²Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6−7 significant figures of precision. Usually white dwarf models without²²Ne are fit to the observed periods and other properties. The rms residuals to the ≃150−400 s low-order g-mode periods are typically in the range ofσ_rms≲ 0.3 s, for a fit precision ofσ_rms/P≲ 0.3%. We find average relative period shifts of ΔP/P≃ ±0.5% for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of ΔP/Pdepending on the specific g-mode, abundance of²²Ne, effective temperature, and the mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when²²Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of²²Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including²²Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model-fitting processes.